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A LIMITED STUDY OP THE APPLICATION OF 
CONCENTRATED LOADS TO FIXED EDGE.AE® DISCONTINUOUS TWO-WAY SLABS 
OF REINFORCED CONCRETE, USING THE CENTRAL DIFFERENCE OPERATOR OF 
THE F3ETTE DIFFERENCE APPROXIMATION TO THE BIHARMONIC EQUATION 

By 

Ronnie Ray Henk 
ABSTRACT 

This paper makes a comparison of resulting design moments for a 
particular size, unifoimly loaded, two-way reinforced concrete .slab 
with' various boundary conditions, using both the Central -Difference 
Operator of the Finite Difference Approximation to^ the Bihamonic 
Equation and Method 3 of ACI Bulletin 5l8-63- The 'atterrpt 'pf this 
study was to determine whether the large grid Finite Difference ‘ 
Analysis would reasonably approximate the results of ACI 31S-63 
Method 3 for a uniformly loaded slab. 

Due to the reasonably close correlation of the results} of the* 
methods the Finite Difference Analysis, with the aid of the Elastic 
Curve Plot, was used' to determine the design moments of, the .slab under 
the influence of a concentrated load. A concentrated load in Finite 
Difference Analysis is defined as being uniformly distributed over an 
eq,uivalent grid area, with the center of the loaded area and the loaded 
grid point being coincident. 



A method of estimating the xE^imom allowable, symmetrically located 
concentrated load is presented, although the results are unproven by 
testing. 

A step-by-step account of the solution to each problem is presented. 
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W. LICT OF SYMBOLS 


A 

A-total 

C 

^subscript 


D 

DL 

h 

LL 

m 

M 


P 

pcf 

psf 

<1 


area of the slab over which the live load is appli 
total slab area 

g.h^ 

an assigned constant for the problem equal to 
moment coefficient for two-way slabs as given 

in tables. 1, 2, and 3 of ACI 3l8-65> Section A2003- 
Coefficients have identifying indexes, such as 


neg’ % DL> • ' • 

flexural rigidity of the plate equal to 


Et5 
12(1 - 


dead load 
grid spacing 
live load 

ratio of short span -to long span for two-way slabs 
moment, in ft-lb. Identifying indexes refer to 
direction, sign, and type of loading, such as 


^ pos LL’ ^ neg DL? * 
an assigned constant for the problem' equal to 



lb per cubic ft 
lb per square ft 

uniform load for entire slab, has identifying indexes 
such as qj3p and qjj^ 

concentrated total live load distributed over one grid 


area, equal to 


JL 

h2 



q. 


t 


R 

t 

T 

w 

W 

x,Y 

X 

Y 


concentrated total live load distributed over 

T 4. W 

areas', equal to 

(3h)2 

an assigned constant for the prohlem equal to 

slab thickness 

allowable concentrated load 

deflection 

total live load 

rectangular coordinate axes 

length of clear ^an in short direction 

length of clear span in long direction 

Poisson's ratio 


nine grid 



D 



V. IHTSODUCTIOM 


Justification for this Study .: ACI Bulletin 318 - 63 ' has "been 

extahlished as the code which "provides minimum requiremen-t;s .for thi 
design and construction of reinforced concrete or composite' structural 
elements of any structure erected' under the requirements of the general 
building code of which this code forms a part." It presents. the require- 
ments for the design of a two-way reinforced concrete slab with varying 
boundary conditions, but only for the application of a uniform load. 

In practice, a uniform loading condition is seldom the usage condition. 
The ACI Bulletin 318-63 and authors of hooks of reinforced concrete 
design occasionally mention that other techniques can be used to analyze 
a two-way slab for non-unifoim loading conditions. They also indicate 
that even though' these methods give approximate answers, the answers 
are often within the realm of reasonably predicting the action of a 
slab structure. 

A design engineer is often faced with the problem of determining 
the capabili'fcy of a two--way slab to withstand a large concentrated load. 
Since ACI 318-63 does not handle -this loading condition,, and since 
ACI 318-63 is. the basis for the reinforced concrete design, one must 
be certain that, any other analytical method used will produce reasonably 
comparable answers. The method used should also be fairly straight- 
forward so that the work may be done efficiently. 

Intent of this Study : To a limited degree, this study will inves- 

tigate the use of the Central Difference Operator of the Finite 
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Difference ApproxinLation -"to' tlie Bihannonic Equation to ascertain "whether 
the resulting' moments are comparable to ACI '^lQ-63 Method 3 * A large 
grid,, ^suitable for de^sk calculator solution, will be used to establish 
,a suitable efficiency in the analysis. 

This study will begin by considering a particular two-way rein- 
forced concrete slab of. dimensions 20 ft ty 20 ft, and determining the 
m ax i m iitn moments for each of three cases of boundary conditions: 

Case "A, " all edges fixed; Case "B, " two opposite edges fixed, two 
opposite edges pinned; Case ”C,^' all edges pinned. The initial deter- 
mination will be made utilizing a uniform loading condition and analy- 
zing by both ACI 518-63 Method 3 and by the Finite Difference Approxi- 
mation to the Biharmonie Equation. A con^arison of the resulting 
maximum moments, will then be made. If a reasonably close correlation 
exists between the results of the two methods, then a reasonable degree 
of accuracy can be expected from the Finite Difference Approximation to 
the Biharmonic Equation to analyze the slab for the application -of -con- 
centrated loads. In the Finite Difference Analysis, a concentrated 
load is defined as being uniformOy distributed over an equivalent grid 
area (h x h) with the center of the loaded area and the loaded grid 
point being coincident. 

In addition to a moment determination for "the application of a 
single concentrated load at the slab center, an atteirpt will be made 
to study the conditions that result from expanding the concentrated 
load until full unifom loading is again achieved. 



Vi. TWO-WAX HIAB. AtlALltSiS - UWIgOKM LOAD 


PesiRn' Criteria 

-(!)■ .The slab is’ to be two-way reinforced concrete of dimensions 


20 ft by 20'ft. 

(2) Minimum slab thickness (t) = greater of ACI 518-63,. 

paragraph 2002(e). (3-1/2 inches or — ^ = O.W5 ft 

= 5.kk inches. Arbitrarily use t = 7 inches to reduce the deflection.) 

(3) Slab dead load = 7/l2 ft by 150 pcf = 87.5 psf 

Live load = 112.5 psf 

Total load = 200.0 psf 

Although ACI 318-63 -Method 3 indicates that the dead load and live load 
are combined in the determination of the negative moment at the fixed 
panel edge, they will be used individually in the negative moment 
calculation for later comparisons. 

(4 ) Eeinforced -concrete- will be considered as having a Poisson' s 
Ratio (p) = 0.15. 

Boundary Conditions 

Case "a" All edges fixed 

Case "B" Two opposite edges fixed - two- opposite edges pinned 

Case "C" All edges pinned 
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Design Asstmiptions 

,(l) Reinforced concrete is elastic, homogeneous, and isotropic. 
Althou^ this is a tremendous oversiiiplification of the exact nature of 
reinforced concrete, to date, few (if of the exact characteris- 

tics can he defined. The reason this difficulty exists is due to the 
vast number of variables that influence its strength. For exainple, 
curing teirperature,. water-cement ratio, strength and gradation of 
aggregate, placement of reinforcing, placement of forms, workmanship, 
etc. , fl.l1 control the strength of reinforced concrete to a certain 
degree. Even though the results obtained will only be an approximation 
of the exact ■ results, this assumption normally gives a reasonable pre- 
diction of the slab action and serves to sin^lify the analysis. 

(2) Plain sections remain plain. 

(3) ' Keutr^ axis is undefofmed. 

(^4-) Deflections are' small in con^arison to the slab thiclmess. 

(•5)' The slab is thin in relation to its linear dimensions. 


Case "A" 


Moment Deterraination. by ACI 3l3-63 Method 3 
All edges fixed 


m = i = 221 = 1.00 

Y 20' 

= Cg • a • x2 

My = • q • 


ACI Case 2 



where 


C = moment coefficient as given in tables 1, 2, and J 
q_ = uniform load 

X and Y = length of respective sides 
Positive Moment Calculation Erom table 2 

q,j3L = D-I- •= 87-5 Psf 


DL - DL ” 

DL " “y DL " (0-01^(87.5)(20)^ = 63O ft-lb 


From table 3 
= L.i;. = 112.3 psf 


IL “ LL ” O-027 


Mx LL = My li = .(0.027)(112.5)(20,)2 = 1215 ft-lb 


Negative Moment CaXculation 


From table 1 
^DL = 

= 112.5 psf 


neg ~ neg “ -O.CA-5 

^ DL neg = \ DL neg = C-0;C*'5)(87.5)(20)2 = -1575 ft-lb 
“x LL neg' = “y.LL neg = (-0.*5)(U2.5) (20)2 ^ .2025 ft-lb 
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SnTrtma.ry j 

% .pos ” 
neg ” 

Case "B" 


Positive 


‘or Case "A” 

% pos = % DL + % LL = 630 ft-lb + 1215 ft-lb = 1845 ft-lb 
% neg = % neg DL + % neg LL = -1575 ft-lb - 2025 ft-lb 

-5600 ft -lb 

Two opposite edges fixed, ACI Case 5 

two opposite edges pinned 

m = _ = 1.0 

= Cg • 4 

' JV= ^ 

Moment Calculation From table 2 

%-L ^ 

•'=A DL = °-°27 

' Cg gg.= 0.018 

Mx DL = '(0.018.)(87.5)(20)^ = 630 ft-lb 
% DL = (0-O2V^CS7.sU.?01^ = 945 ft-lb 

From table 3 
q.gg = 112.5 psf 

^A Lb = 0.032 


Cp li “ 0-027 



Mx ii'= (0.027)(tL12.5)(20)^ = 1215 ft-lb 
My. = ( 0 . 032 )(ll 2 . 5 )(' 20 )^ = IWO ft-lb 



From table 1 

iBL = psf 
•lli = 5 psf 



The ACI Code provides tliat at discontinuous edges- the tending moment 
in the column strips shall be gradually reduced from the full moment 
value of the middle strip to l /5 of these values at the panel edge. 
Therefore, 


% -pos panel edge = 5 ( 18^5 ft-lb) = 615 ft-lb , 

■My poa = My .poa. DL t My p,, = 9^5 ft-lb +'mo ft-lb = 2585 ft-lb 

% neg = % neg DL t My neg. IL ” -2625 ft-lb - 5375 ft-lb = - 6 OOO ft-lb 
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Case "C** All edges pinned ACI Case 1 

m = | = X.O 

= Cg • .q • X2' 

% = • q • 

Positive Moment Calculation From tatle 2 

qni, = 87.5 psf 
q^i = 112-5 psf 

Ca = Cb = 0..036 

Hjj. Bl = ^ DL = (0-036)(87.5 K 20)2 =1260 ft-Vo 

Mx li = Mjr II = (0.036)(U2.5)(20)2 = l620 ft-lb 
Summaiy for Case "c" 

“k pos = “y pos = “y DL + ^ LL = ^260 ft-lb + 1620 ffc-lb = 2880 ft-lb 

Middle strip 

For tlie discontinuous 'edges the moment is l/5 of the middle strip 
moment.' Therefore, 

pos panel edge = ^ pos- panel edge “ j(2880) = 960 ft-lh 
There are' no negative moments involved in this prohl^. 
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Moment Determination ~by Finite Difference Approximation 


^to the Biharmonic Equation 


I’inite Difference*' Approxlanatlon to Case "A” (see Appendix) 


h ,= £ = I 

If 


Boundary Conditions 
All edges fixed 
Therefore 





5,y 5,y| 


due to 

rotation 

restraint 


also 


^ 5 ^ 


^x,0 = ^0,y=%,lf=%y=° 

due to translation 
restraint 

In applying a uniform 
load (q) over the entire 
surface area, conditions 
of symmetry hecome 
•apparent. Due to these 
symmetrical conditions 
small alphabetic letters 
will be assigned to 
points of common deflec- 
tion, to out down on 
writing and make 
symme-try more obvious . 


Superimposing the Pinite Difference Approximation to the Biharmonic 

q h^ 

Equation over each point of differing deflection: Let C = g 
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At grla poiin; <±, d. 

' • 20a, - 8b -^8b - 8b - 8b + 2c + 2c + 2c + 2c = C 
grouping' common terms 

20a - 32b + 8o = C (l) 

At grid point Z,ij' same for l,2j _ 2,3 and 3,2 

I • • i ■ ’ 

. -20b'- 8a-- ;8c -8c + 2b'+2b + b + b = C 

grouping common terms 

-8a + 2 . 6b - l6c = C ’ (2.) 

At grid point 1^1} same for 5^3 

20c - 8b - 8b + 2a + c + c + c + c = C 
grouping common terms 

2a - l6b + 2kc = C (3‘) 

Placing the simultaneous equations in matrix form 


Sum 



( 3 )- 

2 

-l6 

2k 

c , , : iGM-c : 

(Ai 

C} = (2) , 

-8 

26 

- l6 

c : 2+c : 


(!)■ 

20 

-52 

8 

c : - 4 +c ; 


The Cholesky or Grout solution yields (see Appendix) 


Check 


( 3 ) , 2 

-8 

12 

.5000c 

: 5+. 5000c : 

l{t i k} = (2) -8 

-38 [ 

- 2.1053 

-..1516c 

' : -1.1055-.1316G : 

H 

t 

ro 

o 

128 

37.4784 

.2093c 

: 1.0+. 2093c : 
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Back substitution yields 
c = .2093C- 

Td - 2.IO530 = - .1316c 

■b = 2. 1053 (.2093c) - .1316c = .W06C - .,1316c 
b = .3090c 

a - 8b, + 12c = .5c •> 

,a = 8(. 3090c) - 12( .2093c,)-' + .5c = 2.I1-T2OC + .5000c - 2.5116c 
a = .4604C 

Substituting back ^into original equations 

(1) 20('.460te) - 32(. 3090c), ,+ 8{rs093C) should = l.OOC 

9.2080c - 9.8880c + 1.67lAc'= + 10.882lt-C - 9.8880c 

= . 99l|Jt.C ~ l.OOC 

(2) , -8(.4604c) + 26( .,3090c) ,-l6(. 2093c) should = l.OOC 

-3.6832c + 8.0340c - 3.3488c = + 8.0340c - 7.O32OC 

= 1.0020C'*> l.OOC 

(3) 2(.4604C) - 16(. 3090c) + 24(. 2093c) should = 1..00C 

.9208c - 4.9440c + 5.0232c = 5.9440c - 4.9440c 

= l.OOC 


Therefore the solutions satisfy the original equations fairly well. 
■The equations ^or bending moments arel^ 




ax2 . 3y2/ 


My 


. . + u ^ 


W 


hx.‘ 


The Central Finite Difference Approximation for the second partial 






derivative is 
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or in modular form 


© 0 


and 


a^vr 1 I I 

— - » — / -2 Wj^ O modular form 

3y2 h2 ' 


a^ _ 1 


© 

i 

© 


ay® h® 

|Wi,j+l ' 

/ 

^ 

The Maximinn Positive Moment occurs at the center of the slat, point 2^2. 
Therefore superimposing the Central Finite Difference Approximation to 
the second partial derivatives over point 2,2 


_ 1 


== - 2a + 1) = -^(2b - 2a) = -%(.3090 - .4604)C 


ax® h®' 


.1514) = . 

h® h‘^ 


^ if-b - 2a + -b) = ^(2b - 2a) = - 

ay®- h® 


.3028c 


^ pos pos 


- H- ^ *<•»>(- ^)) 


Since 


5028 22(1.15) = .5482 i q ^ = .3482 q.h' 


^DL ! Ili = ^^®-5 psf 


b = 5 = 1 = ^ = 5- 

4 4 4 
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My pos DL = M^ DL *= -5^2(87 .'5 )('25') = 762 ft-lb 

My pos LL = Mj, pos li, .5W2(.112.5.)(25) = 979 ft-lb 

The Maximum' Negative Moments occur at, the center of each fixed panel 
edge, points 0,2; .2,0.;. '4,2 and 2,1)-. 

" ' \ ' • 

.^x neg. 0y2 ~ Jieg 4-.^2 = ^ neg 2>0 ” neg 2,4 ” “ 

^ , 2(0) + b.) = a = 2 il3O|0Cl ^ ,6^ 

ay2 h^ h^ h"^ 

Therefore 

Mueg .6180 q.' = - .6180 q 

h'^ D 

“x neg DL = My neg DL “ - .6l80(87-5)(25). = - 1552 ft-lb 

Mx neg LL = My neg LL “ - -6l80(ll2.5.)(25) = - 1758 ft-lb 

Summary for Case "A" 

Mx neg ~ My neg = My neg DL t My neg LL " 5°90 ft-lb 



Mx pos ” My pos - My pos dL My pos LL 1741 ft-lb 



Finite Difference A-pproxlkation to Case "B" (see Appendix) 


h 


- X = I 
“It. If 


Boundary Conditions 


Two opposite edges 
fixed, two opposite 
edges pinned 

Iherefore 



V allowing 




- “Wj / rotation of 
pinned edge 



Wx, 0 = WQ, y = 4 = WA^ y = 0 

due to translation 
restraint 


In applying a uniform 
load (q) over the entir-e 
surface area, conditions 
of symmetry become 
apparent. Due to these 
symmetrical conditions 
small alphabetic letters 
will be assigned to 
points of common deflec- 
tion, to cut down on 
writing and make 
symmetry more obvious . 


Superimposing the Finite Difference Approximation to the Biharmonic 

Equation over each point of differing deflection: Let C * ^ ^ 

■ - ■ n 



At grid point 2,2 


20a - 8b - 8b - 8c - 8c + 2d + 2d + 2d + 2d = C 
grouping common terms 

20a - l6b - l6c + 8d C 

At grid point 1,2/ same for 5,2 

20b -8a-8d-8d+2c + 2c+b-b = C 
grouping common terms , ’ / 

- 8a + 20b + ifc - l6d = C 

At grid point 2, 5j same for 2,1 

20c - 8a - 8d - ’8d + 2b + 2b + c + c = C 
grouping common terms 

‘f 8a + 4b + 22c - 2j6d = C. 

• At grid point l,5j same for 1,1; 3>1 and 3^3 

■20d - 8b' - 8c '2a + d*+ d + d'- d = C 
grouping 'Common terms 

2a - 8b - 8c + 22d = C 

Placing the simultaneous equations in matrix form 

Sum 


(^) 

2 

-8 

-8 

22 

,c‘ 


8+C 

(2) 

-8 

20 


-16 

G 


C 

(5) 

-8 


22 

-16 

C 


2+C 

(1> 

20 

-l6 

-16 

8 

C 


-4+0 
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rne unoj.esity or uroui: soxutiion yxexas ^see Appenaix; 

Check 


" 2 

-4 

-4 

11 

.5000c 

, ; 4+. 5000c 

-8 

-12 

2.3533 

-6.0000 

-.’4167c 

' ; -2 .6667-. 4167c 

-8 

-28 

55.5524' 

, -1.7350 

-.1205c 

: -.7350-. 1205c 

20 

64 

-85.3312 

23.95,04 

. 3084 c_ 

. : 1.0+. 3084c 


Back substitution yields' 
d = - 5084 c 

c - 1.7350a = - .1205c 

c = i.7550(. 5084c )‘ - .1205c = .5351c - .1205c 
c = .4146c 

b’+_ 2.3533c - 6_..0d = -- .4167c 

■b = -.2.3333(.4 i 46 c) +■ 6..0( .3084c) - .4167c = 1.8504c - 1.3841C 
b = '.46630 

• a -4b - 4c + 'lid •=■ .5000c, 
a = 4(. 4663c) + 4(.4l46c) ll-(. 3084c) + .5OOOC 

a = 4.0236c ■- 3.3924c 

a = .6312c 

Substituting' back into original equations 

(1) 20 (. 6312c) - 16 ,(. 4665c) - I 6 (. 4 l 46 c) + 8(. 5084 c) should = l.OOC 
12 . 6240 C - 7'.46 o 8 c - 6.6336c + 2.4672c = 15.0912c - 14 . 0944 c 

= .9968c ~ l.OOC 

(•2) -8(. 6312c) + 20(. 4665c) + 4(.4l46c) - l6(. 3084c) should = l.OOC 

-5.0496c + 9.3260c + 1.6584c - 4.9344c = 10.9844c - 9 -9840c 


, = 1.0004c ~ l.OOC 
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(3.) -8( .6312c) + 4 (. 4663c.) + 22 (. 4 l 46 c) - 16 (. 3084 c) should = l.OOC 
-5.0496c + 1.8652c + 9.1212c - 4 . 9344 c = 10 . 9864 c - 9.9840C 

= 1.0024c ~ l.OOC 

(4) 2(. 6312c).- 8(. 4663c)' - 8(.4l46c) + 22(. 3084c) should S l.OOC 
1.2624c - 3-T304C - 3.3168c ’+ 6.7848c = 8.0472c - 7.0472c 

= l.OOOOC 

Therefore the' solutions satisfy the original equations fairly well. 
The equiations for hending" moments are (see F.D, Approx, for Case "A") 





8x2/ 


The Maximum Positive Moments occur 'at the center of the slah, point 2^2. 

Therefore^ superimposing the Central Finite Difference Approximation to 
the second partial derivatives over point 2,2 


^ - i(h - 2a + h) = ..^(2b - 2a) = - .6312)C 

8x2 1,2 , 1,2 1,2 


= ' 2 |(- . 1649 ) = ■- 

b- 


.3298c 


& « -IL(c - 2a + 0) = -i-(2c - 2a) = ■^{.klk6 - .6312)0 
8y2 h2 fa2 

= .2166) = - 
h2 
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^ K- ^ 

« - BS(- .5298 - .0650) = .59W 4 ^ 

.59“^8 q 

0495) 
5' 

“x pos LL “ .3948(112.5)(25) = 1110 ft-lt> 

My pos DL ” •4827(8 t.5)(25) = 1056 ft-lb 

My pos LL “ -4827(112.5) (25) = 1358 ft-lb 


“y 


pos 


- D 


..l!iy£,(.15)(^.iffi))=.p(. .4352'-. 


. 1 }. 82 T q_ = . 11-827 1 


Since 


9.DL = ®'7'-5 Psf ' <1 lL = Psf 


h = £ = I = 20' 

4 4 ■ 4 


^ pos dL “ ••5948(87.5)'(25) = 864 ft-lb 


The Maximum Negative Moments occxir at the center of each fixed panel 
edge, points 2,0 and 2,4. 
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^ - i(e - 2(0) + e) = ^ = ^(.lnlf6c) 

=< .8292 4 

My neg - - d|-8292 - .8292 ^ q = - .8292 q 

My neg DL ' -8292(87.5)(25.) = - 0.8l4 ft-lU 
My neg LL .8292(ll2.5)(-25) = - 2J52' ft-llD, 

Sunmiary for Case' "B" 

neg “ ^ neg DL ■*■ neg LL ^ ~ 4l46 ft-lL 
pos .= pos DL + Mic pos IL: ft-lb 

My pos = My pos DL + My pos LL “ SlUlj- ft-lb 
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Finite Difference 'Approximation- to Case "G" (see Appendix)' 


^undary Conditions 


h 


- X _ I 
" '■ -U 


All edges pinned 
Therefore 


-c ‘ -b ' -e 



■"-l,y = -n,y': 
■«'x,-l = -'"x,l 
"■"'5,y| 


due to 
allowing 
rotation of 
pinned edge 


"x,,5 




■"xp- ’”’4,y ° 

due to translation 
restraint 


In applying a uniform 
load (q,) oyer the entire 
surface area^. conditions 
of symmetry become 
apparent. Due to these 
-symmetrical conditions 
small alphabetic letters 
•will be assigned to- 
points of common deflec- 
tion^ to cut down on 
writing and make 
symmetry more obvious . 


Superimposing the Finite Difference Approximation to the Biharmonic 

h^ 

Equation over each point of differing deflection: Let ^ “ ‘I- 

At grid point 2^2 

20a - 8b-8b-8b-8b + 2c + 2c + 2c + 2c-C 


grouping common terms 

20a- - 52b + 8c = C 


(1) 
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At grid poiat 2jlj same for 1,2; 2,5 and 3,2 

20b - 8a - 8c - 8c + 2b + 2b + b - b = C 
grouping common terms 

-8a + 24b - l6o = C 


( 2 ) 


At grid point 1,1; same for 1,5; 5^1 5^5 

20c - 8b - 8b + 2a + c + c -c - c = C 
grouping 'ccramon terms 
2a - 16b + 20c = C 


O). 


Placing the simultaneous equations in matrix form 

Sum 

(5) 

^ r c) = (2-) 

(1) 

'The Cholesky or Grout solution yields (see Appendix) • , 

Check 


2 

-16 

20 

: c, 

' : - 6 +c : 

i . 

-8 

20 

-16 

; c 

' ■*. Q • 

20 

-52 

8 

: c' 

i • 

r -Ih-g : 




2 

1 -8 

10 : 

.5000c 

: 3+. 5000c : 

.-8 

- 4 o|^ 

-1.6000 : 

■-.1250c 

: -.6000-. 1250c ; 

20 

128 

12.8000 : 

.5469c 

: i.ot. 5469c : 


Back substitution yields . 

e = .5469c 

b - 1.6000c = - .1250c 

b = 1.6000(. 5469c) - .1250c = .8750c - .1250c 


b = .7500c 
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a - 8b + IGc = .5000c 

a = 8(.750GC.) - 10(.5lf69C) + .5OOOC = 6.5OOOC - 5.*t690C 
a = 1.0310c 

Substituting back into original equations 

(1) 20(1. 0310c) - 32(.T500C) + 8(. 54690) should = l.OOC 

20.6200c - 24.0000C + 4.37520 = 24.9952c - 24.000GC 

= .9952c ~ l.OOC 

(2) -8(1.03100) + 24(. 7500c) - 16(. 5469c) should = l.OOC 

-8.2480C + 18.0000c - 8.7504c = 18.OOOOC - 16.9984c 

= 1.0016c ~ l.OOC 

(3) 2(1.03100) - 16(.7500C) + 20(.5469C) should = l.OOC 

2.0620c - 12.0000c + 10.9380c = I3.OOOOC - 12.0000c 

= l.OOOOC 

Therefore the solutions satisfy the original equations fairly well. ^ 

The equations for bending moments are (see F.D. Approx, for Case "A") _ 



The -Maximum Positive Moments occiu:' at the center 'of the slab, point 2,2. 
Therefore superimposing the Central Tinite Difference Approximatiqn to 
the second partial derivatives over point 2,2 
2 

i-H i(i3, _ 2a + b) = ^(2b - 2a) = -%(-T500 - 1.0310)C 

8x2 j,2 1,2 j^2 

= Sg(. .2810) =■- .5620 

h h2 
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^ » ^(ID - 2a + b) = i(2b - 2a), = - -5620 •% 

3y2 h,2 h2 

My pos = Mx pos “ - I>(- .5620 ^ + (.15) .5620 


« .5620 2S(i.15) = .6465 4 1 ^ 


.6465 <1 


qp,. = 87.5. psf qpp = 112.5 psf h = f = f = 5' 


My pos DL = Mx pos DL - .6465(8T.5)(25) = l^Ht ft-lb 


My pos LL = Mx pos li “ .6465(112.5,)(25) = I818 ft-lb 


There’ is no negative moment involved in. this case. 


SuTTimary for -Case "C” 


M =0 

neg 


pos “ ^ pos " ^ pos DL pos LL " ^252 ft-lb • 
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Comparison of Results for Uni-fona Loading Condition 


Maximm 

Moment 

Grid . 
Point 

ACI 518-65 ■ 
Method 5, 
ft-lh 

Central Finite 
Difference 
Approximation, 
ft-lh 

Finite 
Difference 
Approximation, 
percent of ACI 

Elastic 

Curve 

Plot 

ft-lh 

Case 

"A" 





Mx 

2,2 ■ . 

+1845 

+1741 

9k’.k 


“y 

2,2 

+1845 

+1741 

94.4 


Mx 

M 

y 

e:Sl 

mi 


-5600 

-5090 

85.8 

-5800 

' ! 

Case 

"B" 





' Mx 

/0,2 

Vlt.,2 

+615 

0 



Mx 

2,2 

+1845 

+1974 

107.0 


. My 

2,2 

+2585 

+24l4 

101.2 


“y 

(2,0 

-6000 

-4l46 

69 

■ ' ^5587 ; 

Case 

"C" 






OJ CVJ 

+960 

0 



“y 

r 2,0 

\z,k 

+960 

0 



■ M^ 

2,2 

+2880 

+5252 

, 112 ..2 


”y 

2,2 

+2880 

+5252 

112.2 



Ta'ble 1. Design Moments for Unifom Loading Condition 










The results compare quite favorably with the exception of the maximum 
fixed edge moments in both Case ”A" and Case . This indicates that 
the grid chosen was too large to closely approximate the deflections that 
influence the maximum negative moment. However, since halving the grid 
intervals approximately quadruples the amount of work necessary to obtain 
the solutions to the simultaneous equations, it would appear that a more 
simple method, of approximating these controlling deflections is in order. 
This can be done by sketching the elastic curve of the slab section 
through the point -of interest, using the calculated deflections. The 
inflection point in the elastic curve is then approximately located and 
the deflection and distance from the point under consideration is scaled 
off the curve. (Note - a thin spring steel wire with straight pins at 
plotted deflection points works wonders in approximating the elastic 
curve) . The second partial derivative is then recalculated using the 
new values obtained and the moment is modified accordingly. 



o o 

^ — ! OJ 

o 

’ II II II 

Is ^ ^ 

From elastic curve 

h' ^ .-172X 
w' « .18C 

^ _i_f2w’) = gl-l8gl = 1£- 

ay^ (h')^ (.172X)^ Y 


Therefore 


ro H 
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with 


Since 


Y = 20', 




9-dl “ • • 


\&y neg DL ' 


p ag-iT ) (87:^)Y^ = 
y2 (256)D 

.0l^79(.87.5)-(■'^00) = 


*^x&y neg- LL ~ 


jj ( 12 . 17 ) . (H2.g)Y^ 

y2 (256)D 


- .0ii75(H2.5.)(40O,) = 



■ = 112.5 pr' ■ ■ 

_ (12.17)C87.5)Y^ ■ 

256 ■ . . 

- l 662 -ft-lh 

(12'.17)(.112-5)Y^ 
' . " 256 

- 2158 ft-ib 


’^&,T neg “ ^ neg DL + neg EL " 3800 ft -lb 







with 


Y * 20*, 




Since 


Idl " • 111 “ ^^^-5 p®^ 

M . D ^ _ .0658M8t.5)x2 

y neg JJL ^ 

== - .0698l^(8T.5)(l^00) = -2444 ft- IL 


% neg LL - D (112.^)]^ = . .06984(112.5)1^ 

256D 

«= - .o6984-(U2.5)(4oo) = - 5145 ft-lb 


“y neg - ^ neg DL + % neg LL 558? ft-lL 



VII. TWO-WAY SLAB ANALYSIS - CONCENTRATED LOADS 


An examination of the preceding Table 1 indicates that the results 
obtained by the Central Finite Difference Approximation to the Biharmouic 
Eq.uation along with the Elastic Curve Plot compares,' well with the results 
'Obtained by the analytical_, emperical of experimental methods used to 
establish Method 5 of ACI 518-65* Therefore^ the Finite Difference 
Approximation can be used to analyze a slab for' the’ application of 
concentrated loads and a reasonable degree of accuracy can be expected. 

This stiidy will now explore the application of concentrated loads 
on the three cases previously examined. It will attempt to investigate 
how the maximimi moments are affected as the load is expanded from' 
application over an area ^ by ^ to full uniform loading. 'The full live 
load (q*X*Y) applied to the slab under the uniform^ loading condition 
will still be used^ only it will be applied over the smaller area. For 
this study, the load will be located such .that the symmetry of the . 
deflections is maintained. The- assumptions previously stated still 
apply. 
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Finite Difference Approximation to Case- "A” 
Concentrated Load Over Grid Point 2,2 (see Appendix) 


h 


* X * 

4 


I 

k 


Boundary Conditions 
All' edges lixed \ 



The total live load previously placed 
on the slah will now he concentrated 
over an area h X h at point 2,2. 


Therefore 


^-1,7 - '^1,7 


%,-l = 

; due to 


i rotation 

II 

< 

restraint 

■^x,5 '^x,5J 



also 




= w,, =0 


due to translation 
restraint 


Since' the deflection 
synmetry conditions are 
retained hy the load 
placement, the- same 
small alphabetic letters 
can be assigned to the 
points of common 
deflection. • 


Previous total live load = 112.5 "paf x 20’ x 20'' = ij-5000 lb .= W 



Superimposing the Finite Difference Approximation to the Biharmonic 

h^ 

Equation over each point of differing deflection; Let R = Qa 
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At grid point 2,2 

20a - 8b-8b-8b-8b + 2c + 2c + 2c + 2c=E 
grouping common tenas 

20a - 52b + 8e = E (l). 

At grid point 2,lj same for 1,2; 2,5 and 5,2 

20b -8a-8c~8c + 2b + 2b+b + b = 0 ■ 

grouping common terms 

-8a + 26b - l6c ■= 0 (2) 

At grid point l,ij same for 1,5} 3,1 aud 3,3 

20c - 8b - 8b + 2a + c + c + c + c = G 
grouping common terms 

2a - l6b + .24c = '0 ( 3) 

Placing the simultaneous equations in matrix fom. 

Sum 

( 5)' 2 -l6 2k : 

^ : c} = (2), -8 26 -16 : 

(1) 20 -32 8 ; 

It should be noted that a similartiy exists between the preceding 
simultaneous eqmtions and those obtained under the uniform losing 
consideration. The only difference is in the value of the load placed 
on the grid points. Therefore, a large portion of the. matrix 



remains the same. 
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Check 


L.'^ : k} = 


Back substituting yields 
c = .0267E 
b - 2.1053c = 0 
b = 2.105'3.( .O26TE) 
b = .O562R 

a - 8b + 12c = 0 

a = 8(..0562R) - 12(.026lR) = .ljk97R - .520te 
-a = .I293R 

Substituting back into original equations 

(1) 20(.1295R) - 52(.0562R) + 8 (.026711) snouia = x.uuk 

2.586OR - I.7987H - .2156E = 2.7996E 7 I..7987R 

= I.O609E 1.00E 

(2) -8(,.1293R) + 26(.0562E) - 16(.026TE) sbould'= O' 

-1.0344E + 1.4615E -..k272E = 1.4615E - 1.46 i6b. 

= - .OOOIR 0 

(3) 2(.1293E) - 16(.0562E) + 24(.0267E) should = 0 

.2586E - .8992 E + .6408e = .8994-E - . 8992E 
= -.0002E ~ 0 

Therefore the solutions satisfy the original equations fairly well. 


2 

-8 

• 12 : 

0 

: 5 : 

-a 

-38 

-2.1055 : 

0 

: -1.1053 : 

20 

120 

37.1|•781^ : 

.0267E 

: i.0f.o267E : 



5 ^ 

The equations for bending moments are (see F.D. Approx, for Case* "A" 
uniform load) 





hy^J 


The' Maximum Positive Moments -occur at the-center of the slab, point 2,2. 
Therefore superimposing the Central Finite Difference Approximation to 
the second partial derivatives over point 2,2 

^ == i(-b - 2a + b) =’i(2b’- 2a) = -%'(.0562 - .1295)R 
h2 h2 h‘^ 


= .0751) - - .11^62 i 


^ J-(-b 

3y® 


- 2a ,+ b) = -■ .ll^62 

•h2 


My' pos LL = Mx pos LL - - •l't^2 J + (.15) (- .lh62 

=» ,1462 ^ (1.15) = -l68l — R 

h^' h2 

« .1681 Srr ^ = .I68IW 

h2 ® 


« .1681(45000) 
=< 7565 ft-lb 
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The Maximum l^egative Moments occur at the center of* each fixed panel 
edge, points 0,2} 2,0} 4,2 and 2,4. Although the Finite Difference 

'Approximation analysis' for the unifonu loading condition indicated the 
grid spacing to 'be too coarse to closely approximate the maximum 
negative moments, the negative moments will he calculated for an order 
of magnitude determination and then checked with the Elastic Curve Plot. 


^ neg 0,2 - •^x neg 4,2 My 2,0 ~ ^ neg 2,4 



^ - 2(0) + h.) = -^(.0562R) 

'h-f- h h^ h*^ 


..1124r 


h 


2 


Mx&y neg LL 


D 


/.1124e\ 

\ h2 -7 


.1142 IL. 


Hh^ 

D 


== - .1142(45000) 


. 114 ^ 


«= - 5139 ft-lh 
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Finite Difference Approximation to Case "A" (see Appendix) 


Concentrated Load Over all Internal Grid Points 


h 


X _ Y 
4 ■ 4 


Boundary Conditions 


All edges fixed 


c b 0 



The total live load previously placed on the 
slab will now be distributed over an area 
5h .X 5hj effectively concentrating a portion 
of the load over each internal grid point. 


Therefore 


11 

> 


= ■"x,l| 

due to 


< rotation 

11 

restraint 



also 



’'x,0='"0,y='"x,4 = %y=° 


due to translation 
restraint 


Since the deflection 
synmetry conditions are 
retained by the load 
placement^ the same 
small alphabetic letters 
can be assigned to the 
points of common 
deflection . 


Previous total live load = 112.5 psf 20' X 20* = i^'5000 lb == ¥ 


at 



Superimposing the Finite Difference Approximation to the Biharmonic Equa- 


tion over each point of differing deflection; Let P 


at 


h^ ■u2 

ii- = 5000 ^ 
D D 
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At grid point 2,2 

20a - 8b - 8b - 8b - 8b + 2c + 2c + 2c + 2c = P ' 
grouping common, terms 
20a - 52b + 8c = P 

At -grid point 2,lj same for l,2j 2,3 and 3j2 

20b - 8a - 8c - 8c + 2b + 2b + b + b = P ' 
grouping common terms 

-8a '+ 26 b - 16c = P 


'( 1 ) 


'■(2y 


At grid point 1,1^ same for ■1^5^ 3^5 

20c - 8b - 8b + 2a + c+-c + c + c = P 
, grouping comon terms 

2a - l6b + 24c = P 


(3-) 


Placing the simultaneous equations in matrix form 


(a •: c} = 


Sum 



-16 

24 

p' 

- : 10+p 

-8 

26 

-16 

p 

: 2+p 

20 

-52 

8 

p 

: -4+p 


Note that the matrix above is identical with that of the original 

uniform loading consideration. Therefore the solutions to the simulta- 
neous equations and the deflection coefficients are identical. The 
difference in the moments will come from the different load applied at 


the grid points. 
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Therefore 

a = .4604 p 
b ■= . 509OP 
0 = .2093P 

The equations for bending moments are (see P.D. Appro?:., for Case "A" , 
uniform load) ■ 


Mx = - D| 






The Maximiim Positive Moments - occur at the center of the slah^ point, 2,2. 


_ , ^028p 

^ ^ ' h2 

My pos LL = M, ^,3 lL » - d(- ^ (-15)(- OOPS ^)) 

~ .5W2 IE = .3I1-82 — (5OGO 

h2 d/ 

= I7IH ft -lb 


The Maximum Negative Moments occur at the center of each fixed panel 
edge, points 0,2; 2,0; 4,2 and 2,4. Althou^ the Finite Difference 

Approximation analysis for the uniform loading condition indicated the 
grid spacing to he' too coarse to. closely approximate the maximum negative 
moments, the negative moments will he calculated for an order of magnitude 
determination and then checked with the Elastic Curve Plot. 



ii-0 


0 

neg -0,2 = ^ neg h,.2 ^ neg '2,0 “ “y neg 2,k ~ ~ ^(5^ ^ 

5SS *6l80P 

MxSy .neg LL « * = ~ -6180 ^(jOOO 

5090;ft-l-b 

The Elastic Curve Plot for -this loading' condition is the same as that, 
previously shown for Case "A" uniform load.. Therefore, the curve data 

i I \ ‘ 

at the inflection point is; 

,±i2Z 

■w* « .*18P 



S2w _ I2.I7P 


Therefore 


with 


^&y neg EL ~ 



- 


^12.17P\ 

. J 



p ^ 9000 h^ _ ■^OOOY^ 

D 16d 


M^&y neg LL - - = - 3805 ft-lL 
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Finite Difference Approximation to Case **B*‘ 
Concentrated Load Over Grid Point 2>2 (see Appendix) 



'The total live load previous]^ placed 
on the slah will now be concentrated 
over an area h X h at point 2,2.. 


Boundary Conditions 

Two opposite edges fixed, 
two opposite edges pinned 


Therefore’ 




due • to 
allowing 
' rotation of 
pinned edge 

due to 

rotation 

restraint 


due to translation 
restraint 


.Since the deflection 
symmetry conditions are 
retained by the load 
placement, the same 
small alphabetic letters 
can be assigned to the 
points of common 
deflection. 


Previous total live load = 112.5 psf X 20' X 20' - ^5000 It = W 


Superimposing the Finite Difference Approximation to the Biharmonlc Equa- 

, ^ „ h^ _ 

tion over each point of differing deflection: Let R - — ^ 
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At grid ■point 2, 2 

20a - 8b - 8b - 8c - 8c + 2d + 2d + 2d + 2d = R 
grouping conunon terms 

20a - l6l3 - l6c + 8a = E 

At grid point' lj2j same for 5,2 

20b - 8a - 8d - 8d + 2c + 2c + b- - b = 0 
grouping common terms 

-8a + 20b + 4c - l6d = 0 

At grid point 2^3; same for 2^d. 

20c - 8a - 8d - 8d + 2b + 2b + c + c = 0 
grouping common' terms 

-8a + 4 Td + 22c - l6d = 0 

At grid point 1^5; same for 1,1; 3»1 3j5 

20d - 8b-8c + 2a+d+d + d— d=0 
grouping common terms 

2a - 8b ^ 8c + 22d == 0 

Placing the simultaneous eq.uations. in matrix form 


(4) 

~ 2 

-8 . 

-a 

22 

:p 


Oi 

( 2 ) 

-8 

20 

4 

-16. 

• 

: 0 - 


'0 

( 3 ) 

-8 

4..-. 

22 

-16 

: 0 


2 

:D' 

20 

> -16 ' 

-16 

8 

: 


. 14 +R 


( 1 ) 


( 2 ) 


( 5 ) 





It should be noted that a similarity exists between the preceding 
simultaneous equations and those obtained under the uniform loading 
consideration. The only difference is in the value of the load placed 
on the grid points. Therefore, a large portion of the L^:I^ matrix 
remains the same. 

Check 


" 2 

' -4' 

' -4 

11 : 

0 


■4 

-8 

-12 . 

' 2.3333 

-6.0000' : 

0 ■ 


-2.6667. 

-8 

-28 

55-332.4 

; -1.7550- ; 

0 


-7350 

20 ; 

64 

-85.3312 

23.9504 

.04 i 8 r‘ 


i.o+.o 4 i 8 r 


Back substituting yields 

d =" .04 i 8 b 
c - 1.7550d = 0 
. c = 1.7550(.04 i 8 e) 
c = .O725R 

b + 2.3333c - e.OOOOd = 0 
b = -2.333(.0725R) + 6.000(.04 i 8 r) 
b = .0816R 

a - 4b - '4c + lid = 0 
a = 4 (.o 8 i 6 r) + 4(.0725 E')' - h(.o 4 i 8 r> 


a = . I566R 



Substituting back into original equations 

(1) 20(.156'6-R) - 16(.08 i 6R) - 16(.0725R), + 8(.01 u 8E) should = l.OOR 

3.1320R - I.3056R - I.1600R + .53!Ar = - 2.1+656K 

■ = I.0008R ~ l.OOR 

(2) -8(.1566 r) + 20(.08 i 6R) + 1)-(.0725 R)-- 16 (. 04 i 8 r) should = 0 

-I.2528R + I.6320R + .2900R - .6688r = 1.9220R - I.9216R - 

= .oo64r ~ 0 

(3) -8(.1566 e) + 4(.08 i 6 r) + 22(.0T25 E) - i 6(.04 i 8 r) should = 0 

-I.2528E + ,3264 r + I.595OE - .6688e = 1.9214R - 1.9216E 

= - .0002R ~ 0 

(4) 2(.1566R) - 8(.08 i 6R) - '8(.0725R) + 22 (.o 4 i 8 r) .should = 0 

.3132R - .6528E - .58OOR+ .9196R = 1.‘2528r - 1.2328R 

= O.OOOOR 

Therefore the solutions satisfy the original equations fairly well. 

The equations for tending moments are (see F. D. Approx, for Case "A" 
uniform load) 


Mx 


= - +-H 

• \ax2 ■ Sy2/ 


•My = 



4 


3x2/ 


The Maximum Positive "Moments bcc,ur at .the center of the slab^ point 2^2. 
Therefore superimposing the Central Finite Difference Approximation to 
the second partial derivatives ovpr point 2j2 ' 



45 


5^ ' ^ 
3x2 


3y2 


-(b - 2a + b) = %2b - .2a) = -%(.08l6 - .1566)E 
.2 h“^ fa2 


= S(- -0750) = - .1500 -| 

h2 IjS 


i(c - 2a + c) = i (2c - 2a> = 4(.0725 - . 1566 .)B 

i2 . 


= ^(.084l) = - .1682 i 
h 


Mk pos LL - - d(- .1500.^'+ (.15)(- .1682) 

. -,(- .I5OO - .0252) = .I752W 

‘ p. 


, »» M752(45000) ^ 

7884 ft-ib 

, y.. ■' . ' • ' 

\'vts LL -1682 J + r.i5)(,-.a5O0)^) 

=» - .(...i682- -■ '.0225') = -1907W 

D 


■ « .1907(45000) 


== 8582 ft -lb 
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The Maximum Negative Moments occur at the center of each fixed panel 
edge, points 2,0 and 2,h-. Although the Finite Difference Approximation 
analysis for the uniform loading condition indicated the grid spacing to 
he too coarse to closely approximate the maximum negative moments, the 
negative moments will he calculated for an order of magnitude determina- 
tion and then checked with the Elastic Curve Plot. 


My neg 



5^w 

^x^ 



2(0) + c)- = % = 4(-°'^25R> 


= .1^^50 \ 


“y 


neg LL 



- .1^50 


D Wh^ 


«= - .ik50(k5000) 







finite Difference Approximation to Case "B“ (see Appendix) 


Concentrated Load Over All Internal -Grid Points 


U 4 



The total live load previously placed on the 
slah will now he' distributed over an area 
3h X effectively concentrating a portion 
of the load over each internal grid point. 


Boundary Conditions 

Two opposite edges fixed 
two opposite edges pinned 

Therefore 

due to 
allowing 
rotation of 
pinned edge 

due to 
rotation 
restraint 




"x,o=%y="x,:4 = %y = ° 

due to translation . 
restraint 


Since the deflection 
symmetry conditions are 
retained by the load 
placement^ the same 
small alphabetic letters 
can be assigned to the 
points of common 
deflection. 


Previous total live load = 112.5 psf X 20* X 20' = 45000 lb = W 


W ^ 45000 = 5000 

(5h)^ 9h^ 


Superimposing the' Finite 'Difference Approximation to' the Blharmpnic Equa- 


tion over each point of differing deflection: 


Let 




^9 


At grid point 2,2 

20a - 8b --Sb - 8c - 8c + 2d + 2d + 2d + 2d = P 
•grouping common terns 

20a - 1613 - l6o +, 8d = ’P 

At grid point lj2j same for 3,2 

20b - 8a - 8d - 8d + 2c + 2c + t - t = P 
grouping common terms 

-8a + 20b-_ + 4c - l6d = P 

At grid point 2^5j same for 2^1 

20c -8a-8d-8d + 2b + 2b + c + c = P 
grouping common terms 

-8a + 4b + 22c - l6d - P 

At grid point 1, same for 1,1; 3^1 saad 5,3 

20d - 8b -8c + 2a + d + d + d- d = P 
grouping common terms 

2a - 8b - 8c -h 22d = P 

Placing the simultaneous equations in matrix form 

Sum 


” 2 

-8 

-8 

22 

P 


8+P 

-8 

20 

4 

-16 

P 


0+P 

-8 

4 

22 

-16 



2+P 

_20 

. -16 

-16 

8 

Pj 


-4+P 


(!■) 


( 2 ) 


) 


) 


(A i 0} = 
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Note that the {a;c} matrix above is ideniiical with that of the original 
uniform loading consideration. Therefore the solutions to the simulta- 
neous equations 'and the deflection coefficients are identical. The 
difference in the moments 'Will* come, from the 'different lo'ad applied at 
the grid points/' 

Therefore 

a = .6512P 
b = .4665P 
c = .ia46P 
d = .5084P 

The equations for bending moments are (see F.D. Approx, for Case "A" 
uniform load) 



The Maxi'TrmTn Positive Moments occur at the center of the slab^ point 2^2 




% pos LL ' 


b(. 


.3298 ^ + (.15)(- 


A552)2g 



=> 1974 ft-lb 
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Mypos LL- - e(- A332.^ + (,15)(- .3298a^j ^ 

- .4827 ^ P = -W27 ^ ( 5000 .^^ 

» 24l4 ft-lb 

The Maximum MesatlTe Moments occur at the center of each fixed panel 
edge, points 2,0 and 2,4. Although the Finite Difference Approximation 
analysis for the uniform loading condition indicated the grid spacing to 
he too coarse to closely approximate- the maximum negative moments, the 
negative moments will he calculated for an order . of magnitude determina- 
tion and then checked with the Elastic Curve .Plpt. 



» - 4l46 ft-lh 

The Elastic Curve Plot for this loading condition is the same as that 
previously shown for Case "B" uniform load. Therefore, the curve data 
, at the inflection point is; 
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h'.»= :iT7y ‘ 
w' <= .23P 

^ 17-88P ■ 

Sy2 y2 

with 

, Y 5000h2 5000 y2 

" 1).’ D “ 16 d 

« - 5588 ft-Ib 
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Finite Difference Approxlmatloii to Case "C” 
Concentrated -Load Over Grid Point 2,2 (see Appendix) 


h 


1 = 1 
"4 1 + 


Boundary Conditions 



The total live load previously placed 
on the slat will now he concentrated 
over an area h X h at point 2,2. 


All edges pinned 
Therefore 


^-i,y 

“■^i,y 




due to 

Wx,-1 

"-Wx,l 

•allowing 



^ rotation of 

""5,y 

= -^5,y| 

pinned edge 


II 



due to translation 
restraint 


Since the deflection 
symmetry, conditions are 
retained hy the load 
placement, the same 
small alphabetic letters 
can be assigned to the 
points of -common 
deflection. 


Previous -total live load = 112.5psf X 20' X 20' - UjOOO lb = W 


W 

‘'"“I? 


Superimposing the Finite 'Difference Approximation to the Biharmonic Equa- 

* ■ • Wh^ 

tion over each point of differing deflection: Let R = qa -^ = 
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At grid point 2,2 

20a -8b'-8b-8b-8b + 2c + 2c + 2c + 2c=R 
grouping common terms 

20a - 32b + 8c = E (l) 

At grid point 2,lj same for 1,2; 2,3 and 

20b-8a-8c-8c + 2b + 2b+b-b = 0 
grouping common terms 

-8a + 24b - 16C = 0 (2) 

At grid point 1,1; same for 3/1 and. 5/3 

20c - 8b-8b + 2a + c + c- c- c = 0 
grouping, common terms 

2a - l6b + 20c = 0 (5) 

Placing the simultaneous eqiiations in matrix form 

Sum 



(3) r 2 

-i6 

20 

0 ; 6 : 

(a : 

CO 

1 

II 

24 

-l6 

0 : 0 ; 


(1) [_20 

-32 

8 ' 

R. : -4+e 


It should be noted that a similarl’ty exists between the preceaxng 
simultaneous eq,uations and .those obtained under the uniform loading 
consideration. The. only difference is -in the Value of the load placed ■ 
on the .grid points.. Therefore, a large portion of matrix 
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raieck 



2 

-8 

10 ' ; 

' 0 

: 3 

: k} = 

'-8 

-40' 

-1.6090 : 

0 

; -.6000 


20 

; 120 

, 12.80 : 

.0781R, 

: 1.0+.0781R 


Back substituting 'yields 
c = .O78IR 
b - 1.6000c f 0 . 
b = 1.6000(.0781R) 

b = .I25OR 

a - 8b + 10c = 0 

■ a = 8(.1250R) - 10(-.07.81R) = l.OOOOK - .7812R' 
a = .2188R 

Substituting back into original equations 

(1) 20(.2188R) - 32(.1250E) + 8'(.0781R) should = l.OOR 

4.376OR - 1^.0000R + .62lt8R = 5.OOO8R - it.OOOOR 

= I.OOO8R ~ l.OOR 

(2) -8(.2188R) + 24(.1250R) - 16(.0781R) should = 0 

-1.T50UR + 3.OOOOR - I.2496R = 3-OOOOE - 3.OOOOE 

= O.OOOOE 

(3 ) 2(.2188e) - 16(..1250R) + 20(.0T81R) should = 0 

.1)-376R - 2.0000R + I.562OE = I.9996R - 2.0000R 

= - .0004R ~ 0 


Iherefore the solutions satisfy the original equations fairly well. 
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The eq.uations for "bending momeni^s -are (see-^F.D. Approx, for Case "A" 
uniform load) 



The Maximum Positive- Moments occur at the center or tne s±ao, poinu 
Therefore superimposing the C^entral Finite Difference Approximation to 
the second partial derivatives over point 2,2 

^ . 2a + b) = ^(2b - 2a) = -§r('1250 - .2188)E 

5x2 ^2 jj2 jj2 


■ = J(- .0958) = - .1876 J 

^ » i(b - 2a + b) = - .1876 4 
Sy h ^ 


. My pos LL = M, pos LL - - .1876 ^ + (.15)(- .1876)^) 


»= .1876 SB(i.i 5) = .2157 ^ E 

h2 h 


= .2157 = .2157(^5000) 

>»*=• D 


« 9706 ft-Ib 


There is no negative moment involved in this case. 
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Finite Difference Approximation 1:0 Case "C” (see Appendix) 
Concentrated Load .Oyer All Internal Grid Points ■ ' 


h 


= X = 
4 


2- 

4 


-0 -b -c 



The total live load preTiously placed on the 
slab will now be distributed over an area 
3h X 3b, effectively concentrating a portion 
of the load over each internal grid point. 


Boundary Conditions 

\ ~ ~ 

All edges pinned 
Therefore 


^-l,y--^l,y'' 
’”x, -1 ~ 1 

^5,y 


due to 
allowing 
rotation of 
pinned edge 


W r- 




due to translation 
restraint 


= 0 


Since the deflection 
symmetry conditions are 
retained hy the load 
placement^ the same 
small alphabetic letters 
can he assigned to the 
points of common 
deflection. 


Previous total live load = 112.5 psf X 20* X 20' - ii-5000 lb = W 


W _ 1^5000 ^ 5000 

(3h)2 9h^ b^ 


Superimposing the Finite Difference Approxin^tion to the Biharmonic Ec^ua- 

,2 

tion over each point of differing deflection: Let P = It = 5000 — 
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At grid point 2^2 

20a - 8b - 8b - '8b - fth + + 2r- + 2a + 2c =.P ’ 

grouping common terms 
20a - 32b + 8c = P 


(. 1 )' 


At grid point -2^1; same for l,2j, 2,5 and. 5,2 

20b - 8a - 8c - 8c + 2b + 2b + b - b = T 
■grouping common terms 

-8a + 24b - l6c = P 

At grid point 1^1; same, for 1,3/ 5/-1 3^^ 3/5 

20c - 8b - 8b + 2a + 'c4-c - c - c — P 
grouping common terms- 
2a - l6b + 20c - P 

Placing the simultaneous equations in matrix form 


Sum 


(3) 

2 

-16 

20 

P 

: ; 

A i e) = (- 2 ) 

-8 

2k 

-16 

P 

: , OfP : 

(1) ^ 

20 

-32 

8 

P' 

: -i(+P‘ : 


( 2 ) 


.(5) 


Hote that the matrix above is identical with that of the original 

uniform loading consideration. Therefore the solutions to the simulta- 
neous equations and the deflection coefficients are identical. Hie 
difference in. the moments will come from the different loads applied at 


the grid points . 
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Therefore 

a = I.O3IOP 
h = .75OOP 
c = .5U69P 

The equations for hend-ing moments are (see F.D. Approx, for Case "A" 
uniform load) 





+ [i 




The Maximum Positive Moments occur at the center of the slah, point 2,-2'. 


^ . .5620 

9x2 1,2 


^ «= - .5620 ■£- 

9y2 h2 


^ pos LL ” “x pos 


li - of- .5620 4 + (•15)(- -5620) 4'' 

\ h h ; 


"« .5620 D (1..15) 


»= .6463 4 (5000) 4 

h2 D 


« 5232 ft-lb 

There is no negative moment involved in this case. 



VIII . SUMMABY OF RESULTS OBTAMED BY THE FINITE DIFgEBEH'CE APPROXIMATIOM 


Case "A" 


A = area 
loaded 


( 5 h )2 

(n )2 

a/a total 

1.0000 

.5625 

.0625 

Positive Moment point 2,2 

%&y pos DL 
^x&y pos LL 
^x&v pos 

+762 

+979 

+i 7 ia 

+762 

+ 17 *H 

+2503 

+762 

+7565 

+8527 

Neeative Moment (center of fixed panel edge) 

^x&y neg DL 
Mx&y neg LL 
neg 

-1662* 

-2138* 

-3800* 

-1662* 

-3803* 

-51+65 

-1662* 

-5139 

-6801 


Case "C" 


A = area 
loaded 

(l+h )2 

( 3 h)^ 

(h)^ 

a/a total 

1.0000 

■ 5625 

.0625 

Positive Moment point 2,2 

l^x&y pos DL 
^&y pos LL 
^&y pos 


+l 4 ll+ 

+3232 

+ 1 + 61+6 

+iini+ 

+9706 

+11120 


Case "B" 


A * area 
loaded 

(l+h)2 

( 3 h)^ 

(h)2 

a/a total 

1.0000 

.5625 

.0625 

Positive Moment point 2,2 

Mx pos DL 
pos LL 
pos 

My pos DL 
Jfy pos LL 
Ify pos 

+ 861 + 

+1110 

+1971+ 

+1056 

+1358 

+ 21 + 11 + 

+864 

+1974 

+2818 

+1056 

+ 24 i 4 

+3470 

+864 

+7884 

+8748 

+1056 

+8582 

+9638 

Negative Moment (center of fixed panel edge) 

^ neg DL 
neg LL 
% neg 


- 2444 * 

-5588* 

-8032 

- 2444 * 

-6525 

-8969 


An * Indicates that the results were obtained using the Elastic Curve Plot 

Table 2 . Maximum Moments Obtained by the Finite Difference Approximation for the 


Three Conditions of Loading 
















!g LI (f-t 


0 .1 

.2 .5 A .5 -6 .7 -8 •? 1.0- 

total 

Figure 2.- Case 

"A" Live Load Mcments Versus Ratio or Loaa 
Area to Total- Area. , 
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A study of the plots of the Live Load Moments versus. A/A total 
leads to an interesting possibility for approximating a maximum concen- 
trated live load that could be placed on the slab without exceeding the 
ACI design moment. 

In the cases examined the design live load is 3mown to be 112.5 psf» 
Therefore the total live load for this particular slab is 
112.5 psf' X 20* X 20' = il-5000 lb. Suppose we are considering placing 
a load at the center of a Case "A" type slab, to bear over an area 
'5 ft X '5 ft. a/a total = (5 ft X 5 ft)/(20 ft x 20’ ft) =. .0625.. In 
Case "A", the ACI design live load moment is 1215 ft-lb. The plot of 
Case "A" shows that a 4-5000' lb live load, at A/A total = .0625, yields 
a moment of 75^5 ft-lb. 

Therefore the approximate allowable concentrated load 

(T) = (-45000 Ib) 

= 7220 lb- 

Since the equation for the moment has been determined in the 
preceding analysis for the load over h X h and A/A total “ - 062 ^, the 
new moment can be checked. 

Case "A" load over h X h Mpog = .I68IT 

= .1681(7220) = 1215 ft-lb 

Therefore the ACI design moment, has not been exceeded but it should 
be noted that this concentrated load does not allow for the placement of 
any additional live load on the slab. 
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After the maximum concentrated load for the positive’ moment has 
been established^ it is necessary to check the maximum allowable load 
related to the negative moment. 

2025 

In a manner similar to above, T ® (^5000 lb) * 17720 lb 

5159 


However, the lesser of the two allowable loads .calcixlated would be 
the maximum load that could be placed at the slab center without ' 
exceeding the- ACI design moment. 

It is important to recall that the load plabed on the slab, for the 
analysis, produced a symmetry in the deflections'. This- symmetry must 
be maintained to use the curves in determining the maximum' allowable 
concentrated load. It would not be maintained if, for example, , the* 
total load was placed one .grid point away from the slab center. ,In 
a case of that nature it would be necessary to use the Finite Difference 
Approximation to analyze the slab with the load in that particular 
position. 

The curve of the Live Load Moments versus A/A total indicates that 
the maximum concentrated load could be determined for all values' of 
A/A total* this study only three points were used to determine the 

curve. Thus, it should be realized that- the curve between these points 
is only roughly approximated. To obtain additional points for this 
curve, it would be necessary to reduce the grid size. ’ A reduction in 
the grid, size' would increase the nimiber of simultaneous -equations 
necessary to solve but it would also increase the accuracy of the 


approximation . 



EC. COMCLUSIOM 


Although the large grid network selected for the analysis tends 
to lessen the accuracy of the approximation, it also effectively 
reduces the number of simultaneous equations which are necessary. to 
solve. A further reduction in the number, of simultaneous equations 
was obtained from symmetry, by identifying grid points of common deflec- 
tion. Table 1 shows that, for a taxiformly loaded slab, the Finite- 
Difference Approximation to the Biharmonie Equation and A&I 518-63 
Method 3 yield congiarable results for the positive moment. Ei the 
area of the fixed panel edges, the Finite Difference Approximation 
did not yield comparable answers. This -was because a grid point did 
not fall near the inflection point of -the elastic curve. It was 
determined that the location of the inflection- point could be -reason- 
ably approximated by plotting the calculated deflections and sketching 
the elastic curve through these points. The negative moment was then 
recalculated using the data from this Elastic Curve Plot. Therefore, 
the, combination of the large grid Finite Difference Approximation to 
the Biharmonie Equation and the Elastic Curve Plot produced an efficient 
and reasonably close approximation to the ACI requirements. 

Since a reasonable comparison -was obtained between the two methods, 
the Finite Difference Approximation to the Biharmonie Equation was- used 
to determine the design moments for the slab under the influence of 
concentrated loads. Table 2 and Figures 2, 3, and 4 show how the 
maximum moments -vary as the total live load is expanded from, application 
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over an area ^ X ^ a ftOX uniform load. In Chapter VIII^ a method 
is shown for approximating the maximum allowable concentrated live load 
which can be applied to the slab without exceeding the ACI design 
moments. It should be noted that for this study the curves in 
figures 2^ 5^ and were drawn through only three points. Therefore^ 
any value picked from the curves, other than those calculated, must be 
considered only as a veiy rough approximation. It should also be noted 
that although the maximum allowable concentrated load does not produce 
calculated moments which exceed the ACI design moments, the results have 
not been proven conclusively by testing. 
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FmilE DIFFERENCE APPROXIMA.TIOI1 TO THE BIHASMOMIC EQUATION 


The equation of equilihrlum for a homogeneous, isotropic flat 
plate in terms of its deflection is. 

D V^w = q(x,y) 


where 

q(x,y) = loading function 
D = Et5 

12(1 - h^) 

E = Young's modulus of elasticity 

t = thickness of plate 

10. = Poisson's ratio 

'^■w = hiharmonic equation 

, a** ■ 

The bihamonxc operator V* = + ^2^2 ^5” 


^ = V^(^) = V^f- 


Laplacian' 

Operator 


The hiharmonic operator operates on the deflection (w) as does, the. 


Laplacian operator. Thus 



Ijow, — 5 is represented "by the Central Finite Difference Approximation 






r±,3-i 
1 j o, 

jy2 ~i,2<- i,j 


Iherefore, V^w is represented 'by 




In modular form 




Now let ^ = V^. Thus, '^'w = = V^j^, 

Since this is the same differential equation as for the deflection, it 


can he concluded that 
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Each of the encircled modules in represents the entire module 
within the 3 brackets of . therefore, |^23 must be superimposed on 
each module of which follows: 
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Therefore, the Central Finite Difference Approximation to the Biharmonic 
Equation is represented by the module 


© 

© © © 

0 0 0 

h^ 

0 0 © 
, © 



With this Finite Difference Approximation module centrally superimposed 
over each point of different deflection on the plate grid system, the 
simultaneous equations are developed in the form 
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13 

Ctiolesky or Grout Method for Solving Sijiultaneous Equations ^ 
Given a set of simultaneous equations in the form 


iOC r C = 0 > 


we -want to get this set of simultaneous equations in a solved fashion 
such that TX - K = 0 where T = upper unit triangular matrix ■ i 


12 


"1 

0 1 
ff 0 

0 .0 
_0 0 

and K is a column vector' 


%5 : 

1 t^ . 


‘'m 


^2n ■ 




O' 0 ... 1 t(„_i)nj 

0 0 ... 0 1 J 


K = 


2 


K . 


To' obtain this solution for the given simultaneous equations^ multipiy 
IX - K = 0 by a lower triangular matrix L such that 




A necessary condition for this operation is that both T and L must 
be' non-singular. Therefore, 


l{0X - K.}.= 0- 
or 

L3X - LK = 0 


Comparing this to original problon, we see that 

LT = A 
LK = C 

Partitioning these matrices to aid in their direct solution 

(a : c } = L (t : k) . 

Expanding these matrices to their coefficients 
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aqi 

ai2 

a-13- • 

• aqji|ci 

agi 

®-22 

agj. 

•®'2n;'°2 

aj]_ 

aj2 

ajj.. 

.aj^j-cj 

%il 

®n2 

3^y: 

•®nn;°n 


^11 

0 

0 

.0 

Zgx 

222 

0 

.0 • 


bs 

253.;. 

;o 

^nl 

^n2 

1^3. . 

• ^nn 


1 tq2i llj-; -lln 

’"I- 

^ ' •'*^2n 

kg 

00' 1 •••tjn, 

kj. 

• • ° ^ ^(n-l)n 

0 0 0 0 1'' 
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The solution of the coefficients of the L and T matrices follows a 
systematic procedure as shown in equation form: 




il 


^13 


^11 


r^ 3 -l 

I 


^±3 ” ®i3 ■ Zj- 

r=l 


'*^i 3 I 


■ii 


"13 


r=i-l 

X 

r=l 




One of the great advantages to using this method for any hand operation 
is that it lends itself to a Check column. 

By summing up all coefficients in each row of the matrix, 

placing each sum ad 3 S'Cent to the row it represents, then a Check column 
equal to all the coefficients in the matrix is obtained by the 

equation 


r=i-l 


Cheeky 


ISuaij^ 


I. 


Checlc^ 


r*l 
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Therefore, as a final matrix form we have 






Sum 






Check 

^■11 

^■12 

^■15 • 



Sumi 


^11 

iti2 t;^5 . 


>^1 


, Check']_ • 

■®21 

®22 

^•23 • 

■•®2n:°2 


Smog 


l2i 

Ig2' ■|‘feg3 


kg 


Checkg 






• 

= . 


I 32 ^331. 

'[ 

^3 


Clieckj ! 

®hl 



• ®nn:°n 


Sumji 


\ 


^nn 

Xn 


, Checkjj 


Note that with this method of comhining the L and T matrix th^ 
1.0's in the diagonal of the T matrix are not written in. Therefore, 
when comparing the results of the Check column, +1.0 must he added to 
the coefficients of the {tIk) matrix to obtain the proper check. 

By back substituting, the solution to the simultaneous equations 
are obtained. That is. 


Xn = 


^(n-1) '’^(n-l)nXn " ’^(n-1) 

^(nv2) ■*=(n-2)(n-l)X(n-l) * ‘'=(n-2)nXn = ^^(n-2) 

etc. 

As a final check, after the con^ilete solution has been obtained, the 
values for Xj_ should be substituted back into all of the original 



